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Abstract 

We present the most recent set of world averages for 
D Q -D° mixing and CP violation parameters, as obtained 
by the Heavy Flavor Averaging Group from a global fit 
to various measurements. The values obtained for the 
mixing parameters when allowing for CP violation are 
x = (0.98 ±8Jt)% and y = (0.83 ± 0.16)%; the signif- 
icance of mixing is 10.2(7. There is no evidence for CP 
violation at the current level of sensitivity. 

Introduction 

In 2006, the Heavy Flavor Averaging Group (HFAG) [ 1 1 
convened a new subgroup to calculate world average (WA) 
values of charm mixing and CP violation (CPV) parame- 
ters Q. Since that time, D°-D mixing has been observed, 
and a wealth of mixing and CPV results have appeared. 
The HFAG charm group has calculated several sets of WA 
values, updating old averages as new results have become 
available. This paper presents the most recent set of aver- 
ages, i.e., those based on results that appeared in preprint 
form by the summer of 2009. 

Mixing in the B° and heavy flavor systems is gov- 
erned by the short-distance box diagram. In the D° system, 
however, this diagram is doubly-Cabibbo-suppressed (rel- 
ative to amplitudes dominating the decay width) and also 
GIM-suppressed. Thus, the short-distance mixing rate is 
tiny, and D°-D mixing is expected to be dominated by 
long-distance processes. These are difficult to calculate re- 
liably, and theoretical estimates for D°-D° mixing range 
over 2-3 orders of magnitude J3JIU- 
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that for q = p, D x is CP-odd and D 2 is CP-even. The 
parameters_P+ = \A{D° -> f)/A{D° -> /)| 2 , R~ = 
\A(D°^f)/A(D°^ f)\ 2 , and0 = Argfe/p). 

To obtain WA values of x,y,S,\q/p\, and 4>, we per- 
form a global fit to 28 measured observables. These 
observables are from measurements of D° — > K + (rv, 

D° -> K + K-/-K+7T-, D° -> K+TT-, D° -> K+TT-TT , 

D° -> K° s tt+tt-, and D° -► K+K-R° decays Q, 
and from double-tagged branching fractions measured in 
e + e~ — > ^(3770) — > DD reactions. To fit these observ- 
ables, we must include an additional strong phase 5 Kn ^ 
(see below). For D° — > K + ir~ decays, we combine 
R + and R~ into parameters R D = (R + + R~)/2 and 
A D = (R+ - R-)/(R+ + R~). Correlations among ob- 
servables are accounted for by using covariance matrices 
provided by the experimental collaborations. 

With the exception of the ^(3770) — > DD measure- 
ments, all methods identify the flavor of the D° or D° 
when produced by reconstructing the decay D* + — > P°7r + 
or D*~ — > D°ir~; the charge of the accompanying pion 
identifies the D flavor. For signal decays, M D , — M D „ — 
M + = Q w 6 MeV, which is relatively close to the 
threshold. Thus, analyses typically require that the recon- 
structed Q be small to suppress backgrounds. For time- 
dependent measurements, the D° decay time is calculated 
/ and 7J° — > / are, respec- as i^/p) x M D0 , where £ is the distance between the D* 

and D° decay vertices and p is the D° momentum. The 
D* vertex position is taken to be either the primary vertex 
position (pp experiments) or else is calculated from the in- 
tersection of the D° momentum vector with the beam-spot 
profile (e + e~ experiments). 
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Input Observables 

The global fit determines central values and errors for 
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All input values are listed in Table [T] The values for ob- 
servables R M = (x 2 + y 2 )/2 0, y CP Q, and A r Q 
are HFAG WA values [8|. They are calculated as weighted 
averages of measurements, taking into account correlations 
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Table 3: Results of the global fit for the cases of no CPV 
0.732 ± 2.890 ± 1.030 % and all-CP V-allowed. 



3.420 ± 1.390 ± 0.740 % 



-1.200 ± 2.500 ± 1.400 % 



-0.500 ± 1.000 ± 0.800 ' 



1.310 ± 0.320 ± 0.250 % 



0.110 ± 0.610 ± 0.520 • 



1.160 ± 0.220 ± 0.180 % 



1.107 ± 0.217 % 
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Figure 1: WA value of y CP as calculated from D° 
K+K~, L 
ments Q. 



K+K~, D° -► vr+Tr-, and P° -► K + K~K S measure- 



among systematic errors and sometimes also statistical er- 
rors. As an example, the weighted average for y CP is 
shown in Fig. Q] The values of observables from D° — > 
K% ir+ir- decays for no-CPV are HFAG WA val- 
ues 1 8 1, but for the CPV^-allowed case only Belle measure- 
ments are available. The D° — > K + ir~ results used [10] 
are from Belle, Babar, and CDF, as these results have much 
greater precision than earlier ones. The D — ► K + ir~ir° re- 
sults are from Babar [ 1 1 1, and the i/>(3770) — > DD results 
are from CLEOc ITU . 

The relationships between the observables and the fitted 
parameters are listed in Table [2] For each set of correlated 
observables, we construct a difference vector V. For exam- 
ple, V = (Ax, Ay, A\q/p\,A<f>) for D° -> K° s tt+tt" de- 
cays, where A represents the difference between the mea- 
sured value and the fitted parameter value. The contri- 
bution of a set of observables to the \ 2 is calculated as 
V ■ ■ V T , where M _1 is the inverse of the covari- 

ance matrix for the measurement. All covariance matrices 
used are listed in TableQ] 

Fit results 

The global fit uses MINUIT with the MIGRAD mini- 
mizer, and all errors are obtained from MINOS. Three sep- 
arate fits are performed: (a) assuming CP conservation 
(A D and (f> are fixed to zero, \q/p\ is fixed to one); (b) 
assuming no direct CPV (A D is fixed to zero); and (c) 
allowing full CPV (all parameters floated). Results from 
the first and last fits are listed in Table [3] For the CPV- 
allowed fit, individual contributions to the x 2 are listed in 
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0.336 ± 0.008 


0.337 ±0.009 
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Table 4: Contributions to the x 2 (CPV-allowed fit) . 



Observable 


x 2 


Ex 2 


VCP 




1.85 


1.85 


A r 




0.15 


2.00 






0.23 


2.23 






2.49 


4.73 






0.00 
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R u {K+l-v) 


0.03 


5.42 


X K+TT-W° 




2.94 


8.36 
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10.04 


R M /y/R D / y 


/ R^'cos5(CLEOc) 


5.72 


15.76 


R+/x' 2 +/y'+ 


(Babar) 


2.74 


18.50 


R-/x' 2 -/y'~ 


(Babar) 


2.01 


20.51 


R+/x' 2 +/y' + 


(Belle) 


3.72 


24.23 


R-/x' 2 -/y'- 


(Belle) 


1.28 


25.51 


R D /x' 2 /v' (CDF) 


0.75 


26.26 



Table g] The total x 2 is 26.3 for 28 - 8 = 20 degrees of 
freedom; this corresponds to a confidence level of 0.16. 

Confidence contours in the two dimensions (x,y) and 
(\q/p\ , <j>) are obtained by letting, for any point in the two- 
dimensional plane, all other fitted parameters take their pre- 
ferred values. The resulting la-5a contours are shown in 
Fig. [2] for the CP-conserving case, and in Fig. [3] for the 
CPy-allowed case. The contours are determined from 
the increase of the x 2 above the minimum value (Xmm)- 
One observes that the (x, y) contours for no-CPT^ and for 
CPy-allowed are almost identical. In the latter case, the 
X 2 at the no-mixing point (x, y) = (0, 0) is 1 10 units above 
the minimum value; this difference corresponds to a confi- 
dence level of 10.2cr. Thus, no mixing is excluded at this 
high level. In the (\q/p\, <j>) plot, the no-CPV point (1,0) 
is within the la contour; thus the data is consistent with 
CP conservation. 

One-dimensional confidence curves for individual pa- 
rameters are obtained by letting, for any value of the pa- 
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Table 1: Input values used for the global fit, from Refs. 161171 l9l HUl HT1 [121 . 



Observable 


Value 


Comment 




VCP 
A T 


(1.107 ±0.217)% 
(0.123 ± 0.248)% 


WA D° — i 


■K+K-fir+iT- and D° -> K+ K~ K° s results (8) 


x (no CPV) 


(0.811 ± 0.334)% 






y (no CPV) 


(0.309 ± 0.281)% 


No CPV: 




\q/p\ (no direct CPV) 


0.95 ± 0.22^0 


WAD°-i 


K° s n+TT- results [8] 


4> (no direct CPV) 


(-0.035 ± 0.19 ± 0.09) rad 







y 
\q/p\ 



(0.81 ±0.30 



+0.13 w 
0.17> /o 

(0.37±0.25t°;^°)% 
0.86 ± 0.30j;°;J° 
(-0.244 ± 0.31 ± 0.09) rad 



CPV-allowed: 
Belle D° 
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results. Con'elation coefficients: 
-0.007 -0.255a 0.216 
1 -0.019a -0.280 

-0.019a 1 -0.128a 

-0.280 -0.128a 1 



Note: a 



-0.007 
-0.255a 
0.216 

(l<z/p| + l) 2 /2isa variable transformation factor 



R M 


(0.0130 ± 0.0269)% 


WAD°^K+£-u results |8 | 


x" 

y" 


(2.61 toil ±0.39)% 
(-0-06 toil ±0.34)% 


Babar D° — > K~*~ir~ n° result. Correlation coefficient = —0.75. 
Note: x" = x cos 5„ ±i/sin<5„ , y" = y cos 5 „ — a:sini5„ 

KlTTT K7T7T ^ K 7T7T K 7T7T 




(0.199 ±0.173 ±0.0)% 
(-5.207 ± 5.571 ± 2.737)% 
(-2.395 ± 1.739 ± 0.938)% 
(8.878 ± 3.369 ± 1.579)% 



CLEOc results from "double-tagged" branching fractions 
measured in -0(3770) — » DD decays. Con'elation coefficients: 
1 -0.0644 0.0072 0.0607 

-0.0644 1 -0.3172 -0.8331 

0.0072 -0.3172 1 0.3893 
0.0607 -0.8331 0.3893 1 
Note: the only external input to these fit results are 
branching fractions. 



R D 

y' + 



(0.303 ± 0.0189)% 
(-0.024 ± 0.052)% 
(0.98 ± 0.78)% 




results. Con'elation coefficients: 
0.77 -0.87 
1 -0.94 
-0.94 1 



A 
x 
V 



D 
12- 



(-2.1 ±5.4)% 
(-0.020 ± 0.050)% 
(0.96 ± 0.75)% 



Babar D° —> K + tt results; conelation coefficients same as above. 



X l2 + 

y' + 



(0.364 ± 0.018)% 
(0.032 ± 0.037)% 
(-0.12 ±0.58)% 



Belle D° — » K + ir results. Conelation coefficients: 

1 0.655 -0.834 

0.655 1 -0.909 

-0.834 -0.909 1 



A 
x 
V 



D 
12- 



(2.3 ± 4.7)% 
(0.006 ± 0.034)% 
(0.20 ± 0.54)% 



Belle D° — » K + tt results; conelation coefficients same as above. 



(0.304 ± 0.055)% 
(-0.012 ±0.035)% 
(0.85 ± 0.76)% 



CDF D° — » A'+7r results. Conelation coefficients: 

1 0.923 -0.971 

0.923 1 -0.984 

-0.971 -0.984 1 
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Table 2: Left: decay modes used to determine the parameters x, y, S, R D , A D , \q/p\, <j>, and S K7rn . Middle: observables 
measured for each decay mode. Right: the relationships between the observables and the fitted parameters. 



Decay Mode 


Observables 


Relationship 


D° -> K+K- /tt+tt- /K+K-K® 


VCP 
A T 


2 y C p = (\i/p\ + \Pll\)y co& 4> - (\q/p\ - \p/q\)xsinij> 
2A r = (\l/p\ ~ \p/<l\)yc°8<t> - (llM + \p/q\)xBincf> 


D°^K% TT+TT- 


X 

y 

Wp\ 

<t> 




D°-*K+i-v 


R M 


R M = (x 2 +y 2 )/2 


D"^K + 7T-TT° 

(Dalitz plot analysis) 


x" 

y" 


x" = x cos 8 tj- + V sin 8 „ 

K 7T7T a K 7T7T 

y" = y cos 8 — x sin 8 „ 

a & Kir it Kir it 


"Double-tagged" branching fractions 
measured in i/j(3770) — * DD decays 


R M 

y 

R D 
^Rq cos 8 


R M = (x 2 +y 2 )/2 




R+, R~ 
x' 2 +, x' 2 ~ 

y'+, y'~ 


R D = {R+ +R~)/2 

A D = (R+ -R-)/(R+ +R-) 

x 1 = x cos 8 + y sin 8 

y' = y cos 8 — x sin 8 

A M = (|,j/ P | 4 - l)/(k/p| 4 + 1) 

x'± = [(1± A M )/(1 T A M )] 1 / 4 (x' cos0±j/ sin<« 

y* = [(l±A M )/(l T A M )] 1 / 4 (y'cos<t> T x> S in<t>) 



Summary 

We summarize the fit results listed in Tableland shown 
in Figs.[3]and|4]as follows: 

• the experimental data consistently indicate that D° 
mesons undergo mixing. The no-mixing point 
(x, y) = (0, 0) is excluded at 10.2a. The parameter x 
differs from zero by 3.2a, and the parameter y differs 
from zero by 4.8a. The effect is presumably domi- 
nated by long-distance processes, which are difficult 
to calculate. 

• Since y CP is positive, the CP-even state is shorter- 
lived, as in the K°-K° system. However, since x is 
also positive, the CP-even state is heavier, unlike in 
the K°-K° system. 

• The strong phase difference 6 is probably not small: 
the fitted value is (26.4 

• There is no evidence yet for CPV in the D°-D sys- 
tem. Observing CPV at the current level of sensitiv- 
ity would indicate new physics. 

Acknowledgments 

We thank the organizers of the CHARM 2009 workshop 
for excellent hospitality and for a stimulating scientific pro- 
gram. 




-1 -0.5 0.5 1 1.5 



x (%) 



Figure 2: Two-dimensional contours for mixing parameters 
(x, y), for no CPV. 



rameter, all other fitted parameters take their preferred val- 
ues. The resulting functions Ax 2 — X 2 — Xmm are snown 
in Fig. [4] The points where A\ 2 — 3.84 determine 95% 
C.L. intervals for the parameters, as shown in the figure. 
These intervals are listed in Table[3] 
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Figure 4: The function A\ 2 = X 2 — Xmin f° r parameters x, y, 6, S K7r7r , \q/p\, and <f>. The points where A% 2 = 3.84 
(denoted by the dashed horizontal line) determine a 95% C.L. interval. 
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